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The main features in iron-pnictide superconductors are summarized as (i) the orthorhombic tran- 
sition accompanied by remarkable softening of shear modulus, (ii) high-Tc superconductivity close to 
the orthorhombic phase, and (iii) stripe-type magnetic order induced by orthorhombicity. To present 
a unified explanation for them, we analyze the multiorbital Hubbard-Holstein model with Fe-ion 
optical phonons based on the orbital fluctuation theory. In the random-phase-approximation (RPA) , 
a small electron-phonon coupling constant (A ~ 0.2) is enough to produce large orbital (^charge 
quadrupole) fluctuations. The most divergent susceptibility is the Orrz-antiferro-quadrupole (AFQ) 
susceptibility, which causes the s-wave superconductivity without sign reversal (s++-wave state). 
At the same time, divergent development of 0^2_y2-ferro-quadrupole (FQ) susceptibility is brought 
by the "two-orbiton process" with respect to the AFQ fluctuations, which is absent in the RPA. 
The derived FQ fluctuations cause the softening of Cee shear modulus, and its long-range-order 
not only triggers the orthorhombic structure transition, but also induces the instability of stripe- 
type antiferro-magnetic state. In other words, the condensation of composite bosons made of two 
orbitons gives rise to the FQ order and structure transition. The theoretically predicted multi- 
orbital-criticality presents a unified explanation for abovementioned features of iron pnictide super- 
conductors. 

PACS numbers: 74.70.Xa, 74.20.-z, 74.20. Rp, 74.25.Kc 



I. INTRODUCTION 

In iron pnictide superconductors^, both spin and or- 
bital degrees of freedom play important roles on various 
electronic properties, such as the high-Tc superconduc- 
tivity, orthorhombic structure transition, and magnetic 
transition. As for the origin of the superconductivity, 
fully-gapped sign-reversing s-wave (s±-wave) state had 
been studied based on spin fluctuation theories^"^. The 
origin of the spin fluctuations is the intra-orbital nest- 
ing and the Coulomb interaction. However, the robust- 
ness of Tc against randomness in iron pnictides indicates 
the absence of sign-reversal in the superconducting (SC) 
gap^"® . 

Later, orbital-fluctuation-mediated s-wave state with- 
out sign reversal (s++-wave) had been proposed based 
on the Hubbard-Holstein (HH) model^"^^. The origin of 
the orbital fluctuations is the inter-orbital nesting and 
the electron-phonon (e-ph) interactions due to non-^i^ 
optical phonons. One of the merits of this scenario is the 
robustness of the s++-wave state against impurities. An- 
other merit is that the close relation between Tc and the 
crystal structure revealed by Lee^^, e.g., Tc becomes the 
highest when the As4 cluster is regular tetrahedron, is au- 
tomatically explained^°. Moreover, orbital- fluctuation- 
mediated s_|--|_-wave state scenario is consistent with the 
large SC gap on the z^-orbital band in Bal22 systems^", 
observed by bulk-sensitive laser ARPES measurement^"^. 

The "resonance-like" hump structure in the neutron 
inelastic scattering-'^'* is frequently explained as the spin- 
resonance due to the sign reversal in the SC gap^^'^^. 



However, experimental hump structure is well repro- 
duced in terms of the s++-wave SC state, rather than 
the s±-wave SC state, by taking the suppression in the 
inelastic scattering 7(w) for |a;| < 3A in the SC state 
(dissipationless mechanism) ^^^8^ distinguish between 
both SC states, measurements of phonon spectral func- 
tion for < 2A would be useful*^. In the normal state, 
prominent non-Fermi liquid transport phenomena in p 
and i?H^° are frequently ascribed to the evidence of spin 
fluctuations^*. However, they are also explained by the 
development of antiferro-orbital fluctuations**. 

In order to identify the mechanism of superconduc- 
tivity, we have to understand the origin of the or- 
dered state, although it is still unsolved in iron pnic- 
tides. For example, in many heavy fermion supercon- 
ductors, the SC phase appears next to the spin-density- 
wave (SDW) state, indicating the occurrence of spin- 
fluctuation- mediated superconductivity. In contrast, the 
ordered state in iron pnictides is not a simple SDW state: 
In fact, the tetragonal to orthorhombic structure tran- 
sition occurs at Tg ^ 100 K, usually above the SDW 
transition temperature Tn^^. In addition, large imbal- 
ance of xz- and j/z-orbitals at the Fermi level and recon- 
struction of the Fermi surfaces (FSs) had been observed 
by angle resolved photoemission spectroscopy (ARPES) 
measurements^'^"^^ . These results indicate that the ferro 
orbital-density-wave (ODW) is the origin of the or- 
thorhombic structure transition. The SDW state may 
originate from the in-plane anisotropy in the exchange in- 
teraction [Jia ^ Jib) associated with the ODW order^^. 

In addition, prominent symmetry breaking C4 —> C2 
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is realized in detwinned 122 systems even above T5 and 
Tn, under very small uniaxial pressure. For example, 
it is recognized as the large in-plane anisotropy in the 
resistivity^^'^^ and the optical conductivity^^ at T* ~ 
200K, which is much higher than T5 and Tn- Moreover, 
the reconstruction of the FSs starts at T* in detwinned 
Ba(Fei_xCox)2As2^^. The discovery of these "electronic 
nematic phase" would indicate that the the ODW fluctu- 
ations dcveolop divergently above T5, and the structure 
transition is (almost) the second-order. 

Recently, prominent softening of shear modulus in un- 
doped and under-doped Ba(Fei_2;Coj;)2As2 had been 
reported by acoustic measurements^°'^^. Yoshizawa et 
al.^^ observed all shear moduli C44, Cge and Ce- Re- 
cently, they had also performed systematic measurement 
for a; = ~ 0.225, and found that only Cge shows 
the prominent softening in both under- and over-doped 
systems^^. Similar results were reported by Goto et al. 
independently^^. This observation can be explained by 
the development of ferro-orbital fluctuations"'^'' or spin- 
nematic fluctuations"^" . Considering large quadrupole- 
strain coupling in iron pnictides, all the observations 
mentioned above suggest the importance of orbital 
physics, and pose a serious challenge for theories of iron 
pnictide superconductors. In fact. Goto et al. have 
shown that Cge is almost independent of the magnetic 
field up to ~ 50T, indicating the non-magnetic origin of 
the softening^^. 

Thus, the main features of the iron-pnictide supercon- 
ductors would be summarized as (i) the orthorhombic (or 
nematic) transition accompanied by remarkable Cqq soft- 
ening, (ii) emergence of high-Tc superconductivity next 
to the orthorhombic phase, and (iii) the stripe- type mag- 
netic order induced by the orthorhombicity. The unified 
explanation has not been achieved as far as we know. 

In this paper, we develop the orbital fluctuation theory 
to explain abovementioned features (i)-(iii) based on the 
random-phase-approximation (RPA) and beyond RPA. 
In the RPA, large O^^z-antiferro-quadrupole (AFQ) fluc- 
tuations are produced by the e-ph interactions, while 
they do not produce the softening of Cee nor C44. If we 
go beyond the RPA, however, we find that the Ox-2_y2- 
ferro-quadrupole (FQ) fluctuations are brought by the 
"two-orbiton process" near the AFQ quantum-critical- 
point (QCP). The induced FQ fluctuations cause the 
softening of Cee, and the commensurate ferro-orbital or- 
der is realized at T = T5. It is predicted that the 
AFQ-QCP is located at the FQ-QCP, which is the end- 
point of the orthorhombic phase. Near this multi-orbital 
QCPs, the superconductivity is mainly caused by the 
AFQ fluctuations, and the orthorhombic transition is 
brought by the FQ fluctuations. Moreover, the stripe- 
type antiferro-magnetic state is induced in the orbital- 
ordered state, since the orbital polarization gives strong 
in-plane anisotropy in the spin-nesting. The present 
study gives a microscopic justification for the anisotropic 
Heisenberg model description in the SDW state^®''^^. 

There is a long history in the study of supercon- 



ductivity due to charge or orbital fluctuations in mul- 
tiorbital systems, starting from the exciton-assisted 
superconductivity^'''"^^. In multiorbital systems, on-site 
Coulomb interaction is composed of the intra-orbital 
term U , inter-orbital one J7', and Hund's or exchange 
term J. Since U is usually larger than U' U — 2 J, spin 
fluctuations induced by U give non s-wave SC states. 
However, Takimoto et a?^° had shown that charge (or 
orbital) fluctuations induced by U' give conventional s- 
wave SC state if the relations U' > U and J ~ are 
assumed. Although this idea was applied to pnictides'*^, 
the relation U' ^ 0.6U is realized in many compounds*^. 
Even if U' ^ 0.6U and J ^ 0.2U, the present authors 
had shown that the orbital-fluctuation-mediated super- 
conductivity is realized by quadrupole-quadrupole inter- 
action mediated by non-Aig phonons, which works as the 
"negative effective exchange J^s for the charge sector."^. 
In iron pnictides, orbital fluctuations develop even when 
the dimensionless e-ph coupling A due to Fe-ion optical 
phonons is just ~ 0.2, according to the RPA-'" and FLEX 
approximation^^. As for As- ion Aig mode**^, orbital fluc- 
tuations develop only when A^ig ~ 1, which is unrealistic 
in iron pnictides. 

Recently, Yanagi et al."^^ added the "orthorhombic 
phonon" to the HH model proposed by the present 
authors^"^^, and studied the ferro-orbital fluctuations. 
However, neither higli-Tc nor structiire transition are 
explained by their theory: First, the "orthorhombic- 
phonon" is acoustic {cOq oc |q|) although it is treated 
as optical in Ref. 43 incorrectly. Their theory belongs 
to the cooperative Jahn- Teller structure transition due 
to acoustic phonons like manganites; see details in Ap- 
pendix A. In this case, the energy-scale of "ferro-orbital 
fluctuations" is too low (~ cOq) to explain high-Tc, since 
experimental Tc is much higher than uiq ~ lOK for 1^1 ~ 
O.Itt^^'^^ SmaU orthorhombicity (a - b)/{a + b)^ 0.003 
in iron pnictides is also inconsistent with the cooperative 
Jahn- Teller scenerio. Second, the derived orbital order is 
"incommensurate;"^'', which is inconsistent with the or- 
thorhombic structure transition and the Cee softening. 

These problems are resolved in the present theory 
since both high-energy AFQ and low-energy "commensu- 
rate" FQ fluctuations develop at the same time, without 
the necessity of fine-tuning model parameters. The for- 
mer (latter) fluctuations give the superconductivity (or- 
thorhombic transition) . 



II. MODEL HAMILTONIAN 

First, we shortly explain the relation between flve- 
orbital and ten-orbital models for iron pnictides. In this 
study, we set x and y axes parallel to the nearest Fe-Fe 
bonds, and represent the z^, .xz, yz, xy, and .t^ — 
d-orbitals in the ccyz-coordinate as 1,2,3,4, and 5, re- 
spectively. The FSs are mainly composed of t2g orbitals 
(1 = 2^4:), although eg orbitals also play non- negligible 
roles in producing orbital fluctuations^". Figure 1 shows 
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the crystal structure of FeAs-plane. Since the As-A (As- 
B) ions form the upper- (lower-) plane, the unit cell con- 
tains Fe-A and Fe-B, which we call the two- iron unit 
cell. Since each Fe-ion contains five-orbitals, the origi- 
nal tight-binding model for iron-pnictides is given as the 
"ten-orbital model" . 

In Refs. 2,42, the authors introduced the gauge trans- 
formation |2,3) — )■ — 12,3) only for Fe-B sites, which 
we call the "unfold-gauge transformation". Due to this 
gauge transformation, the unit cell of the kinetic term is 
halved to become the single-iron unit cell, shown in Fig- 
ure 1. In the obtained "five-orbital model", the kinetic 
term is given as 

C4i^Cj,™<,, (1) 

where i,j denotes the unit cell, l,m = 1 ^ 5 represent 
the d-orbital, and a = ±1 is the spin index, cj^j^ is the 

creation operator of the d-electron, and t^^^ with i ^ j 
{i — j) is the hopping integral (local potential). This five- 
orbital model is convenient to study the Eliashberg gap 
equation^'^. In studying the orbital physics, however, we 
have to keep the fact in mind that the sign of quadrupole 
operators Oxz and Oyz at Fc-B sites are reversed by the 
unfold-gauge transformation. By taking care of this fact, 
we study the softening of shear moduli based on the five- 
orbital model hereafter. In Appendix B, we calculate the 
orbital fluctuations using the original ten-orbital model, 
and make comparison between results of two models. 



two-iron single-iron 
unit cell unit cell 




FIG. 1: (Color online) Crystal structure of FeAs-layer, in 
which the unit cell is given by the two-iron unit cell composed 
of Fe-A, Fe-B, As-A, and As-B. The single-iron unit cell is 
realized by applying the "unfold-gauge transformation" . 

Figure 2 shows the FSs in the (a) five-orbital model 
and the (b) ten-orbital model. The FSs in (a) coincide 
with the FSs in (b) if we fold the former FSs into the two- 
iron Brillouin zone (BZ). We use the hopping parameters 
for LaOFeAs given in Ref. 2. The colors correspond to 
2 (green), 3 (red), and 4 (blue), respectively. The inter- 
orbital nesting between the orbital 2 on FS and orbital 
4 on FS /32 causes the most divergent AFQ fluctuations. 

Next, we introduce the e-ph interaction due to Fe-ion 
Einstein optical modes. The Hamiltonian given in Eq. 



(4) of Refs. 10 is simply rewritten as the following bilin- 
ear form in the a;?/z-coordinate: 

i^e-ph =llY. (0;.< + OL^; + Oly<) , (2) 
i 

where rj = 60e^a^/7-\/3i?pj,_^j,; is the radius of d- 
orbital {e.g., Shannon crystal radius of Fe^"*" is 0.77 A), 

o 

and i?Fc-As ~ 2.4 A. it* is the displacement vector of 
the i-th Fe-ion, and Op (F = xz,yz,xy) is the charge 
quadrupole operator given as 

di^Y.o'r^i^. (3) 

Ini 

where ra] = c]^ ifjCi^ma, and the coefficient is defined 
as o^'™ = 7{l\xz\m) for F = xz, where x — x/r and so 
on. The non-zero coefficients are given as 

o^''=o^'' = V3oi-^ = l, (4) 
-o^f = o^^ = VSol;^ = 1, (5) 
ol^^ = -V3olyy2 = 1. (6) 

Be careful not to confuse Oxz with the a;z-orbital opera- 
tor. Other two quadrupole operators are 0^2 and 0^2 , 

whose coefficients are respectively defined as o^'"^^2 = 
{7/2){l\{x^ - y^)\m) and o'^^ = (7/2^3) (Z|(3i2 - l)\m). 
(They are written as O2 and 0| in literatures.) The 
non-zero coefficients are given as 

c>x2_y2 = -o^^_y2 = -iV3/2)olf_y2 = 1, (7) 
ol'2 = 2o^2^ = 2o^2^ = '~'^z^^ ^ ^ 2/\/3. (8) 

Expect for T — z'^, all the matrix elements of 6r with 
respect to the i2g-orbital (2 ^ 4) are ±1. 

Here, we derived the e-ph interaction based on the 
point-charge model. Although e-ph interaction is also 
induced by the change in the d-p hopping, as discussed 
in Ref. 47, we expect it is small since the weight of 
p-electron on the Fermi surface is just ~ 5% in iron pnic- 
tides. Fortunately, because of the Wigner-Eckart theo- 
rem, the matrix elements of the local quadrupole-phonon 
interaction is always given by the quadrupole operator 
Op , independently of the details of the interaction. Since 
the magnitude of the hexadecapole-phonon interaction is 
{ad/ Rfc-FcY ^ 0-1 times that of the quadrupole-phonon 
interaction, we can safely use Eq. (2). 

Equation (2) means that the displacement produces 
the quadrupole potential Oyz, which causes the scatter- 
ing of electrons between orbitals 2 and 4. The e-ph inter- 
actions in Eq. (2) within t2g-orbitals are shown in Fig. 

3. Then, the phonon- mediated el-el interaction Vj,'j''^^[ is 
obtained by taking the contraction of m*, which gives 
the local phonon Green function D{t) = (rT-uJ^(r)u^(O)) 
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FIG. 2: (Color online) (a) FSs for n = 6.0 in the unfolded 
model. The colors correspond to 2 (green), 3 (red), and 4 
(blue), respectively. The inter-orbital nesting between FS 02 
(green) and FS /?2 (blue) causes the AFQ fluctuations, (b) 
FSs for the original ten-orbital model. 




FIG. 3: (Color online) Inter-orbital scattering processes due 
to the e-ph interaction by {fi = x, y, z) within the t2g- 
orbitals (/ = 2 ~ 4). 



III. RANDOM-PHASE-APPROXIMATION 

Now, we explain the RPA for the five-orbital HH 
model'*^. The irreducible susceptibility in the five orbital 
model is given by 



Xw .mm 



k 



-q) Gl,, {k) 



ill) 



Hi] 



^ is the d electron Green 
= (g,w/), k = (fe,e„), and 



where &{k) = [ie„ -f /i 
function in the orbital basis, q 
e„ = (2n -I- 1)'kT is the fermion Matsubara frequency. 
ji is the chemical potential, and Hi is the kinetic term. 
Then, the susceptibilities for spin and charge sectors in 
the RPA are given as^° 



r iq) = 



X" {q) 



1 - r-x° (q) 
X'{q) 



1 - r-(c.Ox° (q) 

where the bare four-point vertices T"''^ are 



(12) 
(13) 



(/i — x,y,z). By taking the Fourier transformation, we 
obtain 



Diuji) 



(9) 



where uji — 2'kTI is the Boson Matsubara frequency, 
LUo is the optical phonon frequency, and -y/ (u^)o — 
^l/2MFeW£) is the uncertainty in position for Fe ions; 

^/{v?)o = 0.044 A for wd = 0.02 eV^^. Then, V^^^^ is 
expressed as the following quadrupole-quadrupole inter- 
action: 
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J, h = h h = h 
J', h — h ^ h ~ I3 
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15), 



(Ktel 

Here, 



(14) 



(15) 



(16) 



we neglect the 



lIq), and g is the phonon- 



where g{uJi) = g ■ ujl^/iujf 
mediated el-el interaction at zero frequency; g 
in the present point-charge model^^^". 



0.34eV 



/ \ \-^xz,yz,yx 

ladder-diagram for phonon- mediated interaction because 
^ 'oi the relation ujd ^ Whand^ ■ 

In the RPA, the enhancement of the spin susceptibil- 
ity is mainly caused by the intra-orbital Coulomb in- 
teraction [/, using the "intra-orbital nesting". On the 
other hand, the enhancement of x'^ in the present model 
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is caused by the phonon-induced quadrupole-quadrupole 
interaction in Eq. (10), utilizing the "inter-orbital nest- 
ing" in the present model. The SDW (ODW) state is re- 
alized when the spin (charge) Stoner factor Q!s(f,), which is 

the maximum eigenvalue of r*(^^x°(<7, 0), is unity. When 
n = 6.05, the critical value of U is Uc ^ 1.26 eV, and the 
critical value oi g (at U = 0) is gc — 0.233 eV. Smallness 
of gc in iron pnictides originates from the better inter- 
orbital nesting. Hereafter, we set the unit of energy as 
eV unless otherwise noted. 

Here, we introduce the diagonal charge quadrupole sus- 
ceptibilities in the five-orbital model as 

x?iq) - EE«rx«'w(g)or'', (17) 

W nun' 

for r — xz,yz,xy. Their momentum dependence at zero 
frequency is shown in Fig. 4 for ac = 0.98. More gener- 
ally, the quadrupole susceptibility is defined as 

X?,r'i'l) = E E 4'xu'n.n.'{q)oT'"''- (18) 

W mm' 

However, its off-diagonal terms with T ^ T' are negligi- 
bly small in the present model in the "xyz-coordinate" . 
In this approximation, in the absence of Coulomb inter- 
action, the quadrupole susceptibility in the RPA is given 
as 

X?{q)^X?'"{q)/{l-2gx?-'iq)), (19) 

where Xr '' irreducible quadrupole susceptibility. 

Considering the fact that Xw mm' i^) takes large value 
ioT I — m and I' — m' , we obtain Xxfil) ~ ■^x'is 25(9) + 
2X^4,34(9) + (2/3)x;2,i2(9)- Since x^fiq) « 2.5^at q = 
(tt, 0), the critical value of g is gc ^ 0.2 in the present 
model. We stress that the relations Xrr'il) — Xr (9)'^r,r' 
and Eq. (19) holds exactly for q = (0, cj;). We utilize this 
relation in calculating the shear modulus. 

As for the contributions by t2g-orbitals (/ 2 ^ 4), 
X^M « X34,34(9)> X^Aq) « X24,24(9)' and x%i<l) « 
X23,23(9)- In Fig- 4 (a), x^Aq) has the highest peak 
at <7 = (7r,0), which is given by the inter-orbital nesting 
between orbital 3 on FS a2 and orbital 4 on FS (3i in 
the five-orbital model in Fig. 2 (a). Also, XyziQ) has 
the highest peak at q = (0, tt) in Fig. 4 (b), due to the 
inter-orbital nesting between orbital 2 on FS ck2 and the 
orbital 4 on FS f32- We will see that x'^z{yz)(^') modified 
by the unfolding procedure. 

Also, Xxyio) in Fig. 4 (c) is given by the inter-orbital 
nesting between orbital 2 and 3, due to the out-of-plane 
oscillations of Fe-ions. The inter-band and intra-band 
scattering processes produce the enhancement of Xxyio) 
at q = (tt, 0), (0, tt) and q = 0, respectively. We note 
that Xxyio) is not affected by the unfolding procedure. 

Figure 5 shows both Xxzio) at q = (7r,0) and Xxy{^) 
as a function of Uc for n = 6.0 and n = 6.05 given by the 
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FIG. 4: (Color online) Quadrupole susceptibilities in the five- 
orbital model for (a) Xxz(q), (b) XHz(q), and (c) Xxy{<l), 
respectively. The used model parameters are n — 6.05, 
T = 0.05, and g = 0.22 (a^ = 0.98). The correlation length 
in (a) or (b) is derived as ^ = n/Aq ~ 3, where Aq is the 
half-width of the peak. Therefore, we obtain the relations 
~ (1 - ac)-^ and c^^ - 2(1 - ac)'^ - 12^^ 

RPA. We see that Xxzi^) develops divergently in propor- 
tion to (1 — Qfc)^^ oc {gc — g)^^, while xSy(O) shows an 
enhanced but saturated value even at g = gc- 

In Appendix B, we will calculate Xril) in the ten- 
orbital model, and make comparison to Fig. 4 in the 
five-orbital model: Although both results coincide for 
r = xUjZ^jX^ — , they are different for F = xz,yz. 
The reason is that the signs of Oxzjyz at Fe-B sites 
are changed by applying the "unfold-gauge transforma- 
tion". As we will explain in Appendix B, the devel- 
opment of X^^iyJS^i^) in Fig. 4 is the artifact of the 
unfold-gauge transformation. For this reason, the correct 
AFQ susceptibility in the ten-orbital model is given as 
X^./„.(g) = X^;/;r'"'"'(g + (TT,^)). The optical modes 
that give the enhancements of Xyzil) at q = (7r,0) and 
(tt, tt) in the ten-orbital model (q = (0,7r) and (0,0) in 
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FIG. 5: (Color online) Quadrupole susceptibilities as a func- 
tion of Oc for n — 6.0 and n = 6.05 given by RPA. Model 
parameters are U = 0.8 and T = 0.05. We can recognize the 
relation XxziQ) oc (1 — ac)~^ oc (^c — g)~^, which diverges at 
Oc = 1 or g = Qc. 



the five-orbital model) are caused by the in-plane Ux os- 
cillations shown in Fig. 6 (a). Also, the enhancements of 
Xxyio) at <7 = (tTjO) and (0,0) are caused by the out-of- 
plane Uz oscillations in Fig. 6 (b). 
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FIG. 6: (Color online) (a) Fe-ion in-plane optical phonons 
with momentum q — {n,0) and q — (n,n). (b) Fe-ion out- 
of-plane optical phonons with momentum q = {n,0) and q — 
(0,0). 



IV. ACOUSTIC PHONONS 

In previous sections, we studied the e-ph interaction 
due to optical phonons, and calculated the quadrupole 
susceptibilities by the RPA. To obtain the shear mod- 
ulus, we also need the knowledge on the e-ph interac- 
tion due to acoustic phonons with momentum q « 0. In 
fact, shear modulus is proportional to the square of the 
acoustic phonon velocity, which is renormalized by the 
electron-acoustic phonon interaction in the presence of 



strong quadrupole fluctuations. In this section, we derive 
the e-ph interaction due to acoustic phonons with q « 0. 
Hereafter, we use the unit h — 1, and take the nearest- 
neighbor Fe-Fe distance apc-Pc as the unit of length. 

Figure 7 shows the transverse acoustic modes that are 
related to (a) C44, (b) Ce = (Cn - Ci2)/2, and (c) 
Cee^^. Now, we calculate the e-ph interaction based on 
the point-charge model, by following the procedure in 
Ref. 9 for the optical phonons. The quadrupole poten- 
tial energies at Fe-site caused by the transverse acoustic 
phonons in Fig. 7 are given as 



F44 
Ve 



R 



3e2 8 
,4 7^''^^ ' 

Fo-As V o 

3e2 8 



-^Fc-As \/3 



xy ■ UE, 



3e' 



-^Fc-As 



V6(x^ - y^) • {(66, 



(20) 
(21) 
(22) 



for both Fe-A and Fe-B sites, where (x, y, z) is the coor- 
dinates of d-electron. = — u-pc (0 = 44, 66,i?) is 
the relative displacements of the nearest As ions from 
the center Fe ion; (wpe) is the displacement vec- 
tor of the As- (Fe-) ion we are considering from the 
original position. Note that the shear strain tensors 
are given as e44(_E) = U44(£;)/(aFo-Fc/2) = 2u/^^e) and 

£66 = U66/(aFo-Fe/v^) = ^/2uqq. 

The corresponding operators in the ten-orbital model 
are respectively given as 



F44 
Ve 

1>66 



O 



xz ■ W44, 



3e^a 



2„2 



^Fc-As 7%/3 



Oxy ■ Ue, 



Se^al 2^/6 



-"■Fc-As 



7 



Ox2-y2 ■ U66. 



(23) 
(24) 
(25) 



Therefore, the acoustic modes in Fig. 7 (a)-(c) couple 
with the quadrupole susceptibilities at q « 0; Xkz(0)j 
X%{0), and x%-y2{0) for (j) ^ U E, and 66 in the ten- 
orbital model, respectively. 

To study the softening in the five-orbital models, we 
have to perform the "unfold-gauge transformation" for 
Eqs. (23)- (25). Under the gauge transformation, Eqs. 
(24) and (25) are invariant, while Eq. (23) is changed to 



VL 



-"■Fc- 



As 



7%/3 



Oxz ■ U44, 



(26) 



where the — (-I-) sign corresponds to Fe-A (Fe-B) site. 
In the "five-orbital model", therefore, the softening of 
Ce and C66 are caused by Xxy{0,0) and x3^_j^2 (0, 0), 
respectively, while the softening of C44 is caused by 
X^z((7r, tt), 0). Therefore, the softening in shear modu- 
lus (6*66 and C44) does not occur within the RPA^*. 

Next, we derive the effective el-el interaction due to 
fe — > transverse acoustic modes. In the case of fc = 



7 



(a) 



As . 



pi'" 



FeiB 
< 



//k 




FIG. 7; (Color online) Displacement vectors uas and upo in 
the transverse acoustic modes that couple with (a) C44, (b) 
Ce, and (c) Cee- Cee mode corresponds to the orthorhombic 
structure transition. 



k- {1, l)/"\/2 and k <^ 1 shown in Fig. 7 (c), the displace- 
ment operator for the As site at is 



E 



1 



2NMujk 



OfcC 



(27) 



where Ofc and aj[ satisfy the commutation relation 

[flfc: Ofc'] = <5fc,fc'i is the mass of As-ion, and utk = i^fe|fc|; 
Vk is the bare acoustic phonon velocity. The Fourier 

transformation of Ug is given as Uk — \J 2Mujk ^'^^ ~^^'^-k)- 
Then, the Fourier transformation of the acoustic-phonon 
Green function -Dfc(T) = {Ttu\{t)u].{Q)) is 



-Dfc(w„) 



("fc>0, 



(28) 



where (m^)o = {l/2Mujk). 

As understood in Fig. 7 (c), the relative displacement 
with the origin at the Fe-ion, u^, is given as \{us' — Us), 
where Rs' = (0.5, 0.5) and Rs = (-0.5, -0.5) with the 
origin at the Fe-ion. Considering that {Rg' — Rs) ■ k — 
\/2fc, its Fourier transformation is given as 



Uk 



E 



Use 



ik 
V2 



Uk- 



(29) 



£)fc(0)fe2 = 2(l/2MLjfc)feVwfc 

Hee = —gee 



1/Mv1, it is given as 



kk' Al' mm' .aa' 



mm 



^^Ika^l'kcrcl^k'^ 

1 



gee 



7?2 



where B 



3e^ 



Mvl' 



(30) 
(31) 



2V3 



; B = 0.95 eV for 



ad = 0.77 A (^Shannon crystal radius of Fe^+) and 
i?Fo-As = 2.4 A. Therefore, 

566 = ^66Q6!o' (32) 

where Cee,o = is the bare shear modulus, and rjee = 
BRp^_j^^ is the quadrupole-strain coupling constant. If 

we put Vk ^ 0.024 eV A {vk ~ 0.018 eV A) according to 
the first principle study^°, we obtain 366 = 0.1 2eV {gee = 
0.21eV). On the other hand, we obtain g = 0.34eV for 
the Fe-ion optical phonons with = 0.02eV^°. Thus, 
gee/g = 1/2 ^ 1/3 in the present point charge model. 

In the same way, we also derive the el-el interactions 
due to the acoustic phonon with k = k ■ {0, 1), shown in 
Fig. 7 (b). For this mode, the relative displacement Uk 

with Rs (0.5, 0.5) and 
— Rs) ■ k — k, its Fourier 



R^ = (0.5,-0.5). Since {Rs 
transformation is given by 



Uk 



Yuse-^'^^^ 



ik 

jUk. 



(33) 



Then, the phonon-mediated el-el interactions are given 
by the second-order terms of Eq. (24). As a result, the 
el-el interactions due to phonons in Fig. 7 (b) is given as 



Hp 



E 



-gE 

kk' ,11' mm' .( 
^'^Ika'^l'kaC^rnk'. 



W mm' 
^xy xy 



(34) 



where 

3e^ 



B' 



12 



gE 



p2 

^Fe-As 



and B' 



-As 



a-d 



Rfc 



-As 



7V3' 



In the same way, we ob- 



tain 544 = gs- Therefore, g^ = ^0^*0.0 = 0-44?766 

for (j) — 44, E. In conclusion, gs — .944 = OAAgee if Vk is 
equivalent for all modes. 



To calculate the shear modulus, we need the quadrupole 
susceptibility in the "fc-limit" , in which we put uj — 
first, and take the limit fc — ^ later For this pur- 
pose, we derive the effective el-el interaction due to the 
transverse acoustic phonon in the fc-limit. The el-el in- 
teraction due to the phonon in Fig. 7 (c) is given by 
the second-order term of Eq. (25). Using the relation 



V. SOFTENING OF SHEAR MODULI 

A. Softening due to one-orbiton process; the RPA 

Here, we calculate the shear modulus given by the one- 
orbiton process using the RPA. For this purpose, we in- 
troduce the following shear modulus susceptibilities in 
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the five-orbital model, in the absence of e-ph interaction 
due to q « acoustic phonon: 

XE = 2x?,(0,0), (35) 
X44 = 2x?,((^,7r),0), (36) 
X66 = 2x^._^.(0,0), (37) 

where the factor 2 comes from the spin degeneracy. They 
are schematically depicted in Fig. 8 (a). Note that XAi — 
2x^^(0,0) in the ten-orbital model. According to Sec. 2 
in Ref. 51, the shear modulus is given by the second 
derivative of the Free energy with respect to the shear 
strain tensor: The expression for the shear modulus 
((/) = i;,44,66) is^i'S^ 

C4, = C^^o - rilx4>, (38) 

where C^^ = v^p is the bare shear modulus, where 
is the bare acoustic phonon velocity and p is the mass 
density, ry^ is the quadrupole-strain coupling constant 
due to the "acoustic phonon" given in Sec. IV. In Eq. 

(38) , the condition for the structure transition, = 0, 
is satisfied when = (?0^(3> 1). That is, the structure 
transition occurs prior to the divergence of Xit>- 

We can rewrite the expression for given in Eq. (38) 
as follows: 

= C^fll^+g^x^l (39) 

X0 = X0/(1 -5<^X0), (40) 

where = fy^C'^J is the effective el-el interaction due 
to acoustic phonon given in the previous section. In Eq. 

(39) , the condition = corresponds to the divergence 
of X0, since X0 is the total susceptibility including the 
e-ph interactions due to acoustic phonons. 

If we put t/ = for simplicity, Eq. (40) is expressed as 

X44 = X44/(l- (.9 + .944)X44), (41) 

XE = X%/{1-{9 + 9e)x%). (42) 
X66 = X66/(l - ffeeXee)^ (43) 

where the suffix represents the bare susceptibility. Ac- 
cording to X44 in Eq. (41), g in Eq. (19) for Xxz is 
replaced with g + 344 when both optical and acoustic 
phonons are taken into account. Therefore, we have 
to reduce g to g — 544 to keep the charge Stoner fac- 
tor ttc and X^z(Q) unchanged. Considering the relation 
544 ^ gE that we derived in the previous section, we ob- 
tain xe ~ Xe- Then, we conclude that (i) C44 ~ C44^o 
since X44 is seldom enhanced in the RPA. Also, (ii) Ce 
softens to some extent since xS/(0) is weakly enhanced 
as shown in Fig. 5, although the relation Ce — will 
not be satisfied because of the relation xe ^ Xe- 

As for Cee given in Eq.(43), Xee in the present model is 
^ 2eV~^, while we estimate gee — 0.1 ~ 0.2eV. There- 
fore, we expect (iii) Cee Ceefi in the RPA, which is 
inconsistent with experimentally observed large soften- 
ing in Cee"^^- In the RPA, the softening in shear modulus 



(6*66, C44 and Ce) is small according to Eqs. (35)-(37) 
and Figs. 4 (a)-(c). In the next subsection, we analyze 
X66 by taking account of the two-orbiton process that is 
not included in the RPA. 

'^44/66/E 

(a) #y'\/V\/^y#■ 

Or Or 



D(k+q) 




FIG. 8: (Color online) (a) Diagrammatic expression for the 
shear modulus susceptibilities (0 = 44, 66, E) in the RPA. 
(b) The second-order term with respect to the third-order 
anharmonic phonon-phonon interaction A(3)«^U66. This di- 
agram represents the virtual process in which an acoustic 
phonon with q — breaks into two optical phonons con- 
serving the total momentum, (c) Two-orbiton term for the 
shear modulus susceptibilities X6?i which is the irreducible 
susceptibility of Xa;2_y2 (0, 0). This term gives the softening 
of Cee- (d) Dominant contribution for the three-point vertex 
A^2_y2 (r, F; q) for T — xz and yz. (e) Self-energy correction 
due to Xxz{<l)- (f) Second-order Maki-Thompson-type vertex 
correction with respect to Xxzil)- 



B. Softening due to two-orbiton process; the 
Aslamazov-Larkin-type diagram 

In the previous subsection, we studied the softening 
of shear moduli within the RPA. However, the obtained 
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softening is very small since only the AFQ fluctuations 
develop in the present model. In this subsection, we 
analyze xee by taking account of the "two-orbiton pro- 
cess" that is not included in the RPA. In usual cases, 
this higher order process is negligible. However, it gives 
divergent increase in xee, and creates the FQ-QCP near 
the AFQ-QCP. For this reason, the orthorhombic struc- 
ture transition (Cee = 0) is induced by the two-orbiton 
process. 

Before calculating the two-orbiton process, let us con- 
sider the third order anharmonic phonon-phonon cou- 
pling ^(3)u|w66) where Ux is the displacement of Fe- 
lon optical mode in Eq. (2), and uqq is the displace- 
ment of As- ion acoustic node in Eq. (22) or (25). Fig- 
ure 8 (b) shows the second-order-term with respect to 
A(3), where D{q) is the optical phonon Green function in 
Eq. (9), and the factor 1/2 is introduced to cancel the 
overcounting with respect to the upside-down diagrams. 
This term gives a self-energy correction for the acoustic 
phonon Green function. This two phonon process would 
be measurable in Raman spectroscopy. By considering 
the e-ph interaction, each D{q) in (b) is replaced with 
D{q) + {rjD{q))^.X%{q). 

Here, we study the ferro-quadrupole susceptibility due 
to the "two-orbiton term" given by Fig. 8 (c), and an- 
alyze for the softening of Cge- Since the coefficient of 
the anharmonic phonon-phonon coupling A(3) would be 
small in iron pnictides, we instead consider the three- 
point vertex given by electron Green functions in Fig. 
8 (d). The two-orbiton term in Fig. 8 (c) is similar 
to the Aslamazov-Larkin (AL) term for the excess con- 
ductivities due to superconducting fluctuations; such as 
longitudinal^^ and Hall''''* conductivities, and the Nernst 
coefflcient^^. Figure 8 (c) is a virtual process in which one 
Oa;2_j^2-type orbiton with q = breaks into two O^z-type 
orbitons with zero total momentum (±Q). After taking 
the momentum summation, the two-orbiton term (c) in 
2D systems would be strongly enhanced for q = since 
Xyzi^) ^ large peak at finite momentum k. The 
mathematical expression for the two-orbiton process in 
Fig. 8 (c) for X66) which we call x™' given as 



given as 



Ax2_y2{T,r';q) = -2r^Tr{Gfc(e„)6^2_^2Gfe(e„) 
x6rGfe4-q(e„)6r'| , (45) 



where the factor 2 in front of Eq. (45) accounts for the di- 
agrams with reversing three Green functions Gk{f-n) ^■ 
G_fc(— e„). Near the QCP g < gc, the most divergent 
quadrupole susceptibility is Xxz{yz)- Therefore, the dom- 
inant contribution for x™ in Eq. (44) will be given 
by the term with F = F' = xz or yz. After the an- 
alytic continuation, the functional form of X?, 
for q Ri Qxz = (tt, 0) or q w Qy^ = (0, tt) would be 
approximately given as 



X?(q,w + i(5) 



ce 



, (46) 



for F = xz or where ^ is the correlation length and 
Wo is the characteristic energy of the fluctuation. The 
relation a ljq^ holds in the RPA. 

Next, we consider the temperature dependence of ^. 
In the FLEX approximation^* or SCR theory^^, the bare 
susceptibility x° is approximately suppressed as x° — oiT 
{a > 0) due to the thermal fluctuations, which are 
described as the self-energy and Maki-Thompson ver- 
tex corrections. In this case, we obtain Xt'v(QjO) oc 
(1 - 5xL(Q.O) + gaT)-^ cx (T - Tafq)-i "based on 
the RPA, where Tafq = -(1 - gxlziQ,^))/ is the 
transition temperature to the AFQ ordered state. Since 
XxziQ^^) we assume the following relations 



e = 1{T-Tafq)-\ 
coo = I' {T -Tafq), 



(47) 
(48) 



TO 
X66 



^i2gfT J2 Ax2_y2iT,r';q,cJi)Ax2_y2{T,T': 

q,i.,r,r' 

x(l + 2gx?{q,coi)){l + 2gx^,{q,coi)) 

l{2gfT J2 Ax2_y2{T,T';q,uJi)Ax2_y2{T,T': 



q,i,rx' 



xXriq^^i)Xr'iq^^i)^ 



where Ax2^y2(r,T';q,uji) is the three point vertex with 
respect to Vge in Eq- (25) and quadrupole operators Or 
and Or': shown in Fig. 8 (d). We used the relation 



Xr(9) = Xri-l)- When U = 0, A^2_^2 for = is 



where 1,1' are constants. Note that ojq^,^ is temperature 
independent^^. By carrier doping, Tafq changes from 
positive to negative, while other model parameters (c, I, 
and I') would be insensitive to doping. As shown in Fig. 
5, X%{Q,0) ~ 2.4 X (1 - ac)-'^ ~ 12^^^ In the case (i) 
Tafq > 0, the relation luq < T is satisfied near Tafq- 
In the opposite case (ii) Tafq < 0, the relation ojq > 
q,u!^ will hold for wide range of temperatures. Note that 
the present phenomenological model in Eqs. (46)-(48) 
is reproduced by the microscopic calculation by FLEX 
approximation**. As for the spin propagator in cuprate 
q,a;p)ipei'conductors, the relation wq > T (wq < T) holds in 
over-doped (under-doped) systems. 
, , Here, we comment on the self-energy correction and 
^ tiie Maki-Thompson- type vertex correction for xJe', 
shown in Figs. 8 (e) and (f) respectively. The former 
term is included in the FLEX approximation, and it gives 
the Curie- Weiss behavior of XxziQ^ 0) given by Eqs. (46)- 
(48), as reported in Ref. 11 or Ref. 56. The latter term 
would be negligible since its temperature dependence is 
smaller than that of the former term. For this reason. 
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wc concentrate on the two-orbiton term in Fig. 8 (c) 
hereafter. 

From now on, we perform the numerical calculation 
of the two-orbiton process in the case (i), in which the 
relations ^ » 1 and Wq «C T" are realized near the orbital- 
ordered state. In this case, the dominant contribution in 
Eq. (44) comes form the terms with T = V = xz and yz. 
Also, we can safely apply the classical approximation, in 
which the terms with cj; 7^ arc dropped in Eq. (44). 
Under these approximations, Eq. (44) is simplified as 



(25)4t Y,{A,2_y2 {xz, xz; q)x% {q, 0)}' (49) 



To calculate A^2_y2 (T, T';q), we introduce a uniform FQ 

potential term H' = h J2i 0J,2_y2, where h is an infinites- 
imal small constant. Then, the three point vertex is given 
as the following Ward identity''^: 



X^P,(g,0;/i)-x^P,(q,0;0)l (50) 



Ax2_y2{r,r';q) = - 

where Xpp, (<7,a;;;/i) is the "irreducible" quadrupole sus- 
ceptibility with respect to g. In the numerical cal- 
culation, we have to fix fi against the change in h. 
Equation (50) gives the correct three-point vertex even 
for [/ 7^ 0. In the case oi U = 0, Xr r' simply 
given as x^x'i^^^i'^ ^) = -"^Efc,™ Tr{6rG'(fc + q, e„ 
uji;h)dr'G{k,en;h)}. The ohtained A^2_y2{xz,xz;q) for 
J7 = is presented in Fig. 9 (a). Similar result is ob- 
tained for U = 0.8eV. 

According to the functional form of Xsz(q,0) in Eq. 
(46), xle'^/T cx E^{x%{q,0)V « ex (T - Tafq)-^ 
in two-dimensional systems. The numerical result for 
Xee^/T given in Eq. (49) is shown in Fig. 9 (b). The 
obtained result follows the relation xJq' /T ~ 0.1(1 — 
ttc)^^- Since 6^^ ~ {1 — ac)~^, the relation x^q /T ~ 
0.6^^ is verified numerically. 

In the same way, two-orbiton processes for other two 
shear modulus susceptibilities, x™ and x™j are propor- 
tional to the square of At{xz,xz; q) for T = xz and xy, 
respectively. However, they are four orders of magnitude 
smaller than x™ , as recognized in TABLE I: This selec- 
tion rule for Ay is understood as follows: According to 
the relation Tv{OtOxzOxz} = for F = xz,xy, we rec- 
ognize that Axz/xy originates from the off-diagonal terms 
of the Green function G;.„i (l 7^ m.) that is much smaller 
than the diagonal terms. For this reason, the two-orbiton 
process is negligible except for x™- 

Here, we discuss the softening of Cge by taking the two- 
orbiton process into account: According to Eqs. (39) and 
(40), we obtain 

-.-1 



<^66^ = C'eeyi +566X66], 

X66 = Q66 + X™ 

l-566(a66 + X66°)' 



(51) 

(52) 



where agg = 2x6e(0,0). Now, we consider the case (i) 
Tafq > and wq «C T. As we have obtained the relation 



r 




2 2 

X -y 


xz 


yz 


xy 


At {xz, xz, Q) 




-0.60 


1.9 X 10"^ 


-1.0 X 10"^ 


3.2 X 10~* 



TABLE I: Three point vertex Ar{xz, xz, Q) for V = x"^ - y"^, 
xz, yz, and xy. Q = (tt, 0) corresponds to the peak posi- 
tion of X?z('7iO) in the fivc-orbital model. Wc recognize that 
At{xz, xz, Q) ~ 0(1) only for V = x^ —y^', This selection rule 
means that xIae ^ 1- 



Xl^ cx re, we put X66° = bmTl{T - Tafq). Since the 
temperature dependence of oee is small, we obtain 



X66 



a66 + &66 T — (aee/ (a66 + ^66))Iafq 



Ts = Tafq J 



1 - 566 ((166 + '>66) 

1 — 566 «66 



566(066 + ^'66) 



T-Ts 
(> Tafq). 



(53) 
(54) 



Then, the difference between Ts and Tafq, which is 
conventionally denoted as Bjx"''^''''^, is given by Ej^ = 
rs(.g66^66)/(l - 566066) > 0. According to Eq. (53), Eq. 
(51) is rewritten as 

T -Ts 

C66 = (^66,0(1 - 566(o66 + he))-^ — > (55) 

1—1 AFQ 

Here, g^Q = 0.1 ~ 0.2eV and age ~ 2eV~^. We stress 
that Eqs. (53)- (55) are vahd only for loq <C T. 

Finally, we calculate the two-orbiton process analyti- 
cally for general value of ujq/T. Since we cannot apply 
the classical approximation that was used to derive Eq. 
(49), we have to perform the analytic continuation^'' of 
Eq. (44). The obtained expression including the quan- 
tum fiuctuation contribution is given as 

—n{x) 

x2Imxg{q,x + iS)Rex%{q, x + i6), (56) 

where n{x) = {e^^ — 1)~^ is the Bose distribution func- 
tion, and wc put Ax2^y2 outside of the qt-summation since 
its momentum dependence is much smaller than that of 
Xxz- Fii'st, we perform the a;-integration using the fol- 
lowing equations: 



oo 

—=y 



2x 



1 

e-^ — I ^ (2n7r)2 + x"^ ' x ' 

7-2 

dx 



(a2 +a;2)(62 -^a;2)2 26(a 6)2 



(57) 
(58) 



where a,b > 0. Then, the expression after the x- 
integration in Eq. (56) is given as 



c'C' I^ujIT Y.{B,coo + 2nnT)-^ + , (59) 
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(a) (0,71) 




where Bg = 1 + ^^(q — Q)^. Next, we take the q- 
summation ~ 57 /o under the assumption ^ 
1. Then, the q-summation of the second term in Eq. (59) 
is easily obtained as c^^^T/S. Also, the q-summation of 
the first term is given as 



Stt 



^0 

2-itT 



^0 



1 - 



^0 

2tiT 



(60) 



0.03 



FIG. 9: (Color online) (a) Obtained A^-2_y-2 {xz, xz; q) for U — 
0. We put n = 6.05 and T = 0.05. (b) x™/^ given by Eq. 
(49) for [7 = and U = 0.8 as function of Oc. Using the 



relation 6^^ ~ (1 — ac 
obtain X66°/r ~ 0.1(1 ' 



given in the caption of Fig. 4, we 



etc)' 



0.6^ 



where ipix) is di-Gamma function, and the cutoff rimax = 
(1 + ^^7r^)ajo/27rr originates from the fact that the q- 
summation is limited to the region |q| < tt in periodic 
systems. As a result, the final expression for the two- 
orbiton term is 



r 



(61) 



where X 



(2.g) 



An 



{Aj.2_y2(^xz,xz;Q)}'^ . Here, we ver- 
ify Eq. (61) in the opposite two limits: In the case (i) 
luq <C T, the di-Gamma functions in Eq. (61) are negli- 
gible. By applying the relation = 1/{T — Tafq), we 
obtain 



T - Tafq 



(62) 



where h — XL The first line in Eq. (62) coincides with 
Eq. (49) since Eq{x?^(q, 0)}^ = c^£,'^/4:Tt. In the oppo- 
site case (ii) ujq S> T, the term T in the curly brackets in 
Eq. (61) is negligible. Taking the relations ip{x) sa log(x) 
for a; 3> 1 and woC^ oc into account, we obtain 



« Xecoo\og{2 + n^e) 

066 log 



T-T, 



AFQ 



(63) 



where h' — X^^wq. Therefore, in the case (ii) Tafq < 
and Wo ^ T, xgg in Eq. (51) is given by replacing x™ 
with &^glog(7r2|/(T - Tafq)) in Eq. (52). In this case, 
the temperature dependence of x™ is much moderate. 
In Sec. VI B, we will discuss the temperature dependence 
of C66 based on Eq. (61). 

In the above derivation, we have neglected the effect 
of mass-enhancement factor brought by the third point 
vertex. If we take this effect into account, both Eqs. (62) 
and (63) are multiplied by the factor (m*/m)^ = 2^ ^ 3^, 
as we will discuss in Sec. VI B. 



VI. DISCUSSIONS 

A. Why 0^2 y2-FQ fluctuations are the most 
divergent for Tafq > 0? 

In this paper, we have studied the development of 
quadrupole susceptibilities in iron pnictides based on the 
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RPA and beyond the RPA. The main fluctuations in the 
present study is the O^z/yz'-^^^ fluctuations, Xxz/yz^^^' 
which are produced by Fe-ion in-plane optical phonons. 
The acoustic phonons q ^ Q with finite energy also assist 
in producing the AFQ fluctuations; see Eq. (41). We also 
find that the Ox2-y2-FQ fluctuations, X66; are induced by 
the "two-orbiton process" described by the AL-type dia- 
gram in Fig. 8 (c). Here, the anharmonic thrcc-phonon 
coupling is produced by the three-point vertex in Eq. 
(45). The two-orbiton process is important in iron pnic- 
tides because of the two-dimensionality. 

We discuss on the total susceptibility by including the 
electron-acoustic phonon interaction, Xp{q), by taking 
the two-orbiton process into account. For T = xz and 

1 - (.9 + 9n)Xxz[Q) 
Xx-.= (0) = , ""'fot^TToy (65) 

l-ff66(X66(0)+X66 ) 

where xJe' is proportional to the square of the AFQ cor- 
relation length DC XxziQ)- When t/ge = 0> therefore, 
XxziQ) and X%-y2{0) diverges at the same time in pro- 
portion to As gee increases from zero, only X^2_y2{0) 
is enhanced because of the absence of the two-orbiton 
process; Ar{xz,xz;q),Ar{yz,yz;q) <C 1 for F = x^; or 
yz as shown in TABLE I. For this reason, the relation 
Ts > Tafq is universally satisfied even if a fully self- 
consistent calculation is performed. 

As results, both the AFQ fluctuations (=origin of high- 
Tc) and the FQ fluctuations (=origin of shear modulus 
softening) develop at the same time, and latter fluc- 
tuations overcome the former near the Tg. The or- 
thorhombic phase transition in under-doped compounds 
is brought by the divergence of the two-orbiton process 

X66- 

B. Softening of C'ee- comparison between theory 
and experiment 

Is this subsection, we discuss the softening of Cqq in 
under- and over-doped iron pnictides based on the re- 
sults in Sec. V. For this purpose, we first estimate the 
magnitude of the three-point vertex Ax2_y2{xz,xz,Q) 
based on the Ward identity given in Eq. (50). Since 
X?.,x.(Q,0;0) ~ {2g)-\ \Ax2_y2\ ~ {2g)-^/6E, where 
SE is the bandwidth of the xz/yz band. Since 2g ~ 0.5eV 
and SE ^ 2eV according to the band calculations^, and 
considering the effect of band renormalization due to the 
mass-enhancement factor m* /m (= 2 ~ 3)^*, we expect 
|A^2_j^2| ^ l{m*/'m) [eV~^]. This rough estimation is 
consistent with the numerical result in Fig. 9 (a). 

Now, we discus the under-doped case with Tafq > 0. 
In this case, x™/r « X^"^ shown in Eq. (62). Using 
the relations c ~ 12 and 6^^ ~ (1 — ac)~^ as discussed 



in the caption of Fig 4, we obtain X ^ 0.7 and xJe ^ 
0.12(TO*/m)^(l — ac)~^- This estimation is consistent 
with the numerical result Fig. 9 (b) if we put [m* /m) = 
1. 

We also discuss the optimum or over-doped systems 
without structure transition, in which the relation wq ^ 
T is satisfied. In this case, x™ b' login^l/iT-TAFQ))- 
Since the temperature dependence of x™ is moderate, 
X™ would be comparable or smaller than aee- There- 
fore, in over-doped systems, the softening in Cee would 
be much moderate, showing a deviation from the Curie- 
Weiss type form in Eq .(55). 

Now, we analyze the temperature dependence of xie' 
and Cee/Cge.o by using Eq. (61). We can fix the pref- 
actor {2g)'^c^^^{Ax2_y2}^/ATT = X(^ in front of Eq. 
(61) based on the relation xj^ /T = X^"^ for luq <C T: 
We obtain X ~ 0.6 according to Fig. 9 (b). Here- 
after, wc; put X = 5.4 by multiplying the square of the 
mass-enhancement factor, {m*/m)'^ ~ 9. We also put 
wo = l'{T - Tafq) with I' = 2, and ^ = 1{T - Tafq)"' 
with I = 0.086 [cV], which means that ^ ^ 2 when 
T - Tafq = 250K. Using the obtained xle, we plot 
Cee/Cee.o in Fig. 10 (b) based on Eqs. (51) and (52). 
Here, we set gee = 0.17cV and age = 2.5cV~"'^. In the case 
of Tafq = lOOK, we obtain Ejt « 27K. In the FLEX 
approximation^^, Tafq changes from positive to negative 
by carrier doping, while other parameters (X, Z, 1/ and 
age) are insensitive to the doping. Similarly to Fig. 10 
(b) , wc can fit the recent experimental data by Yoshizawa 
et alP for Ba(Fei_xCoa;)2As2 with a; = 0.1, by 
choosing Tafq while other parameters {X , /, V and agg) 
arc fixed. This fact is a strong evidence for the success of 
orbital fiuctuation theory in iron pnictide superconduc- 
tors. 

In the; present paper, we consider that the origin of 
high-Tc is the AFQ fiuctuations. On the other hand, 
Yanagi et al.^^ claimed that high-Tc originates from the 
FQ fluctuations that give the softening in Cgg: In the 
latter mechanism, a rough estimation of T^ is given as 

Te ~ Wcexp(-(1 + /3A)//3A), (66) 

where ujc is the phonon energy relevant for the orbital 
fluctuations, which is just ~ lOK for |qr| ~ O.Itt. /3 = 
1 + 566X66 = Cqq^/Cqq is the enhancement factor due 
to FQ fluctuations^*. However, Cgg.o/C'ee observed in 
optimally-doped Ba(Fe,Co)2As2 is just 1.2^-"-^-'^: Ap- 
parently, such small enhancement cannot reproduce high- 
Tc superconductivity in iron-pnictides. 

In the present study, in contrast, weak softening in 
optimally-doped sample is ascribed to the change in the 
scaling of x™' not to the weakness of AFQ fluctua- 
tions. In fact, the softening is moderate in the case of 
Tafq = -lOK in Fig. 10, while the AFQ correlation 

« 1000/(T[K] -|- 10) is enough to cause the supercon- 
ductivity at Tc ^ 30K. Therefore, moderate softening 
and high-Tc are compatible in the present study. 
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FIG. 10: (Color online) (a) x™ for Tafq = lOOK, 50K, lOK 
and -lOK. given by Eq. (61). (b) Cee/Cee.o for Tafq = 
lOOK, 50K, lOK and -lOK. given by Eqs. (51) and (52). We 
put X = 5.7, I = 0.086, /' = 2, gee = 0.17, and aee = 2.5. 
Cee = is realized when xi? ~ Qm ~ '^ee, which is 3.38 
in the present parameters. Using these same parameters, we 
can fit the recent experimental data by Yoshizawa et al.^^ for 
Ba(Fei_a;Coa;)2As2 for x = ~ 0.1, just by changing Tafq. 



C. Quadrupole-ordered state in under-doped 
compounds 

Here, we consider the orbital or quadrupole ordered 
state in under-doped compounds. In the mean-field ap- 
proximation for the multiorbital Hubbard model for iron 
pnictides^^'^", stripe- type SDW order occurs for U > Uc, 
and weak orbital polarization (uxz 7^ n-yz) is induced 
as the secondary order when the magnetization is large. 
However, in real materials, orthorhombic transition oc- 
curs in the paramagnetic state, and the SDW-order is 
induced in the orthorhombic phase. To solve this prob- 
lem, we studied the multiorbital HH model beyond the 
mean-field theory, and found that the FQ order occurs in 
the paramagnetic state due to the two-orbiton process. 
Fortunately, this FQ order does produce the experimen- 
tally observed SDW order, as we will explain in the next 
subsection. 

As discussed in Sec. VB, the divergence of xee, 
which is the total FQ susceptibility given by both opti- 
cal and acoustic phonons, causes the orthorhombic struc- 
ture transition when Cee = 0. The Ox2_y2-FQ order 
is realized in the orthorhombic phase. The schematic 
quadrupole order is shown in Fig. 11(a). Since 02;2_j,2 « 



713 according to Eq. (7), the order parameter 
y2 > (< 0) corresponds the orbital polarized state 



n2 

with Uxz > TT-yz {rixz < nyz). Figure 11 (b) shows the 
AFQ order brought by the divergence of x'^ziQ)- ^i" 
though the FQ order in (a) would occur earlier, we expect 
the AFQ order in (b) would coexist with the FQ order 
when the structure transition is the weak first order. In 
fact, the reconstruction of the FSs above Tn in detwinned 
Ba(Fei_a;Co2;)2As2^^ would indicate the presence of the 
AFQ order^i. 



(a) 







FIG. 11: (Color onhne) (a 
divergence of 6*55 



0^2_ 2-FQ order given by the 



the — j/^-orbital operator 
the divergence of Xa?z(Q)- 



Be careful not to confuse 0^2_y2 with 



(b) Oa;2-AFQ order brought by 
(c) The correspondence between 
0^2_y2-quadrupole order (Oi2-quadrupole order) and the d- 
orbital with larger occupation number. 



In Fig. 11 (c), we show the correspondence between 
the quadrupole order and the d-wavefunction with larger 
electron occupancy. In the Oa;2_j,2-type quadrupole 
order, the electrons mainly occupy the state \xz) for 
Ox2-y2 > 0, or the state \yz) for Ox2^y2 < 0. In the 
Oajz-type quadrupole order, the electrons mainly occupy 
the state \xy) + \xz) for Oxz > {\xy) — \xz) for Oxz < 0). 
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Finally, we make comparison between the present 
study and the previous work based on the RPA^'^, which 
claims that the divergence of Eq. (51) is caused by large 
566066 ^ 1 while neglecting x™ • However, the obtained 
02.2_j,2-quadrupole order is "incommensurate"'*^. This 
result highlights the importance of the two-orbiton pro- 
cess X™ in order to produce the "q — 0" orthorhombic 
structure transition. 



D. Stripe magnetic order produced by 0^2_y2-FQ 
order 



In under-doped iron pnictides, the collinear-SDW or- 
der is induced in the orthorhombic phase at Tn, which 
is slightly lower temperature than T5. Various expla- 
nations for the origin of this SDW transition had been 
proposed previously. From a strong-coupling scheme, 
square-lattice Heisenberg model with in-plane anisotropy, 
Jia-Jib-J2 model, had been studied^^. According to the 
neutron scattering on CaFe2As2'^'*, the high-energy spin- 
wave dispersion indicates the relation Jia S> Ju in the 
orthorhombic phase. In this case, experimentally ob- 
served staggered spin order along the x-axis (a-axis) is 
expected to be realized. However, such strong in-plane 
anisotropy Jia ^ Jib is surprising, considering the small 
orthorhombicity (a— 6)/(a-t-6) ^ 0.003. These fact would 
indicate the existence of orbital or quadrupole order in 
the orthorhombic phase. 

In this subsection, we study the origin of SDW state 
based on the weak-coupling approach. Hereafter, we as- 
sume that X-axis corresponds to a-axis (longer lattice 
constant) in the orthorhombic phase. In the previous 
subsection, we explained that the two-orbital process 
induces the 0^2_j^2-FQ order in Fig. 11. Note that 
0.j.2_y2 Ki n2 — according to Eq. (7). The corre- 
sponding mean-field is given as 



(a)x«(A£=0) 



ff' = Ai?^(|2)(2|-|3)(3|), 



(67) 



which raises (lowers) the energy-level of orbital 2 (3) by 
AE. In a similar model, the change in the DOS and 
FSs by the orthorhombic potential AE was studied by 
Chen et al.^^. Here, we study the change in the spin 
susceptibility by AE using the RPA. 

We calculate the total spin susceptibility x^iQ^O) = 
Ei.™X?,,;™,™(q,0) for U = 1.1 and g = 0. Figure 12 (a) 
shows the obtained x"!?; 0) = 0; the correspond- 

ing spin Stoner factor is as — 0.87. When AE is finite, 
the four-hold symmetry in x'*(9iO) disappears quickly. 
Figure 12 (b) shows the change in the spin susceptibil- 
ity, x"{q;^E) ~ x'{q-:0), induced by AE = +0M. We 
see that x'*(9iO) increases by +6.5 at q — (0, tt) while 
decreases by —4.0 at q = (''',0). Therefore, magnetic 
frustration is resolved and stripe-SDW order can be in- 
duced by small AE'. 

Figure 12 (c) shows the AE-Uc phase diagram given 
by the RPA, that is, by the mean-field-approximation. 
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FIG. 12: (Color online) (a) x"iQ;AE) for AE = 0. Used 
model parameters are n — 6.05, U — 1.1, g ~ 0, and T = 0.05 
(b) x°(q; AB) - x°{q;0) for AE = 0.04. Model parameters 
are the same as those in (a), (c) Uc as function of AE given 
by the RPA for m* jra = 1. (d) An2 and Ana as function of 
AE for m* jm = 1. Note that An = An2 - Ana. (e) FSs 
for AE = —0.04 for n = 6.05. We use the same color coding 
as in Fig. 2. (f) FSs for AE = -0.08 for n = 6.05. In both 
cases, the best nesting vector is q = (tt, 0). 



It is noteworthy that JJc quickly decreases in proportion 
to \AE\, because of the degeneracy of orbital 2 and 3. 
When U < Uc, the experimental SDW order with mo- 
mentum q — (tt, 0) is realized by the negative potential 
AE that corresponds to n2 > n^. Figure 12 (d) gives the 
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relation between An = An2 — An^ and AE: If we take 
the band-renormalization effect into account, we obtain 
the relation An = —0.85{m*/m)AE. According to (c) 
and (d), we obtain the reduction in Uc due to the FQ 
order is 

AUa = -1.4|An| = -1.2(m7m)|Ai;|. (68) 

Therefore, only few percent An can induce large change 
in Uc that is linear in \An\. According to recent ARPES 
measurement in detwinned BaFe2As2^^, AE ^ — 0.03eV 
in the orthorhombic phase, which corresponds to An — 
+0.026{m*/m) and AUc = -0.036{m* /m) in the present 
five-orbital model. In this case, the realized SDW order 
is q = (7r,0), which is consistent with famous strip- type 
SDW state in mother compounds^^. 

We can show that the SDW temperature Tn also in- 
creases linearly in \An\ based on the Landau theory. The 
free energy in the present problem would be given as 

^^(An) = F(0) + cAn(m^,^o) - mf,^,^), (69) 

where mg is the AF order with momentum Q, and 
F(0) = a.(T-r«)(m2^ ,,^-fm^„^^))+6(m4^^^^+m4^^^P/2+ 

• • • with a, 6 > 0. Then, we obtain ^ + |cAn|/a. 
The present study shows that c < 0, which seems consis- 
tent with the numerical result in Ref. 60. 

Now, we consider the reason why SDW order is pro- 
duced by AE: Figure 12 (e) and (f) shows the change in 
the the FS structure with AE. We can recognize that the 
intra-orbital (orbital 3) nesting between FS a2 and FS /?i 
becomes better, compared to the case of AE = in Fig. 
2 (a). Therefore, the origin of the "FQ-order-induced 
stripe-SDW" is the "anisotropy in the intra-orbital nest- 
ing" caused by small |Aii^| ~ O.OSeV, which corresponds 
to a small orbital polarization |An| ~ 0.026(m*/m). 
This result is consistent with the very small orthorhom- 
bicity (a — b)/{a + b) < 0.003 in the orthorhombic 
state^^. In the strong-coupling description, the origin 
of the stripe-SDW state is the in-plane anisotropy in the 
exchange interaction {Jia ^ Ju)^^'^"^ brought by two- 
orbiton process. 

If we go beyond the RPA, the SDW state will be further 
stabilized by the reduction in the quasiparticle damping 
7 when the FQ-order is established^""^: In fact, in the 
FLEX approximation^^, x'*(q) is suppressed by 7 due to 
strong orbital fluctuations in the normal state. Since the 
orbital fluctuations is suppressed when the AFQ-order 
sets in, the resultant increment in x^il) would stabilize 
the SDW phase. 

E. Summary 

In the present paper, we have studied the realistic five- 
orbital HH model for iron pnictides. In the RPA, only 
the Ox2-AFQ fluctuations develop as shown previously^, 
and therefore the softening of shear moduli (Cee , C44 and 
Ce) cannot be reproduced. In the present study beyond 
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FIG. 13: (Color online) The phase-diagram for iron-pnictide 
superconductors obtained by the present orbital fluctuation 
theory. Ts is the orthorhombic transition temperature (= FQ 
order temperature), and Tn is SDW transition temperature. 
The fact that two QCPs at Tg = and Tafq = almost 
coincide means that novel "multi orbital QCPs" are realized in 
iron pnictides. The left-hand (right-hand) side of the vertical 
dotted line corresponds to Tafq > (Tafq < 0), in which 
the two-orbiton process is relevant (irrelevant). At Tafq, the 
AFQ-order does not occur since it is prevented by the FQ- 
order at Ts- 



the RPA, we revealed that both Oxz-^FQ, fluctuations 
and Oj.2_y2-FQ fluctuations develop at the same time. 
The former and the latter fluctuations are the origins 
of the s++-wave superconductivity and the orthorhom- 
bic structure transition, respectively. The commensurate 
FQ fluctuations are brought by the two-orbiton process 
in Fig. 8 (c) that is dropped in the RPA. [In the mean- 
field theory, the orbital order due to large gge is always 
"incommensurate"'*^.] Fluctuation-induced softening oc- 
curs only in Cgg out of three shear moduli because of the 
orbital selection rule for the three-point vertex. The ori- 
gin of softening would be interpreted as "virtual anhar- 
monicity of lattice vibrations" that is induced by AFQ 
fluctuations; see Fig. 8 (b). Possible quadrupole orders 
in the ordered state are show in Fig. 11. Using the 
two-orbiton term in Eq. (61), we can fit the recent ex- 
perimental data of Cge in Ba(Fei_a;Co2;)2As2'^^ for wide 
range of doping, only by choosing Tafq while other pa- 
rameters are fixed. This fact is a strong evidence for 
the success of orbital fluctuation theory in iron pnictide 
superconductors. 

In addition, we should stress that the stripe-type 
antiferro-magnetic state is realized in the orbital-ordered 
state, since the small orbital polarization (An < 0.05) 
can cause large in-plane anisotropy in the exchange in- 
teraction (Jia 7^ Jib)- Thus, the present study presents 
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a microscopic justification for tlie anisotropic Heisenberg 
model description for the SDW state^^'^^. 

In Fig. 13, we summarize the phase-diagram of iron- 
pnictides given by the present orbital fluctuation the- 
ory beyond the RPA. We stress that Tafq, which is 
determined experimentally from Cge, is positive in the 
under-doped case (T5 > 0) while it is negative in the 
over-doped recognized from Eq. (54) obtained 

in the classical approximation. Especially, « for 
Tafq = 0, consistently with cxpcrimonts"^^'^"^. This re- 
sult indicats that QCPs for AFQ and FQ orders almost 
coincide at the endpoint of the orthorhombic phase. The 
emergence of "multi orbital QCPs" is favorable to the 
orbital-fluctuation-mediated S4-+-wave SC state^"^^. In 
fact, Tafq is derived from experimentally observed Cqq 
as follows: In under-doped systems with T5 > 0, Tafq is 
given the Weiss temperature of Cqq oc {T — Ts)/{T — 9), 
and Tafq = is indeed positive experimentally. In over- 
doped systems, both T5 and 9 are negative (T5 > 9), 
and Cee starts to deviate from the Curie- Weiss behav- 
ior. These experimental results arc the strong cvcdcnce 
for the realization of the two-orbiton process {xJe' ~ 
T/(T— Tafq)) in iron-pnictides. In contrast, in the coop- 
erative Jahn- Teller scenario due to large cjca by Yanagi et 
al^^, the parameter 9 is always negative; see Appendix A 
in detail. Finally, we note that the two-orbiton term x™ 
in the present study is very similar to the bare nematic 
susceptibility Xo.nem in Ref.'^", which is the two-magnon 
term on different sublattices in our terminology. 

In summary, the present study can explain the su- 
perconductivity, orthorhombic transition, and softening 
of Cee due to FQ and AFQ quantum-criticalities. The 
stripe-SDW order is naturally produced by the "or- 
thorhombicity" of the FQ order. These results are strong 
evidence for the realization of the orbital-fluctuation- 
mediated s_|-+-wavc superconductivity in iron pnictides. 
Finally, we stress that the present study enables us to de- 
rive the important parameters in the orbital fluctuation 
model in Eqs. (46)-(48) from the experimental data of 
shear modulus. 
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Appendix A: Structure transitions in itinerant 

electron systems: Quadrupole-quadrupole 
interaction V.S. Cooperative Jahn- Teller effect 

In this paper, we studied the structure transition 
due to Oa;2_j,2-FQ order in iron-pnictides. The ferro 
quadrupole interaction originates from the two-orbiton 
process with respect to Oaj^-AFQ fluctuations, which are 
induced by optical phonons. The 02,2_j^2-FQ fluctuations 
give the softening of the elastic constant Cqq. 

In this section, we present a general theory for the 
structure transitions in itenerant metals, and consider the 
uniqueness of iron pnictides in the next step. The struc- 
ture transition due to ferro-quadrupole order is classifled 
as the "cooperative Jahn- Teller type" or "quadrupole- 
quadrupole interaction type"^^. The elastic constant 
is given by Eqs. (39) and (40), where x<P is the 
quadrupole susceptibility at q = without acoustic 
phonons. Now, we introduce the quadrupole-quadrupole 
interaction g like in Eq. (10), the origin of which is the 
optical phonons or Coulomb interaction. In the RPA, 



— A/0 



X^/(l - 9X4,), where is the bare quadrupole 



susceptibility. Then, is given as 



Cd,.> 



1 - 



9x4, 



(Al) 



In rare-earth metals in which /-electrons are localized, 
x" is proportional to l/T^^'^"^. However, this replace- 
ment is inappropriate in itinerant metals. In nearly ferro- 
quadrupole itinerant metals, |1 — (5'ac + 5')X0l ^ 1 at zero 
temperature. 

Here, we consider the temperature dependence of 
beyond the RPA. In the FLEX approximation^^ or SCR 
theory^^, x^ is replaced with x^ — olT (a > 0) due to the 
thermal fluctuations, which are mainly described as the 
self-energy. In this case, Eq. (Al) becomes 



C41 



T-Ts 
T-9 



(A2) 



where T5 = -(1 - {g^^ + 5)x°)/(5ac + 9)a and 9 = 
~(1 ^ .9X°)/.9CK- Therefore, C4, shows the Curie- Weiss 
behavior in the RPA by taking the thermal fluctuations 
into account. 

First, we consider the case (i) gac ^ g, in which the 
structure transition is driven by quadrupole-quadrupole 
interaction^^. In this case, Ej^ = T5 — 6' « gac/g^o:, 
which is much smaller than T5. (Note that C4, is not 
soften when g^c = 0.) Then, the lattice distortion 
in the ordered state will be very small (^ 1%), like 
in PrRu4Pi2. In the opposite case (ii) ^ac ^ 9, the 
structure transition is driven by cooperative Jahn- Teller 
effect"''^. In this case, 6 w —1/ga, which takes a large 
negative value since a is small; 9 ~ — 300K. Then, the 
lattice distortion in the ordered state should be very large 
(^ 1%), like manganites. In case (ii), the energy of 
the "ferro-orbital fluctuations" induced by the acoustic 
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phonon is very low, and the Tc of the orbital-fluctuation 
superconductivity is still lower. Therefore, emergence of 
high-Tc superconductivity is unlikely in case (ii). 

In the case of Cge in iron pnictides, Oa;2_j^2-0^2_j,2 
interaction g due to optical phonons is absent. Moreover, 
0.j.2_y2-F(^ fluctuations due to Coulomb interaction do 
not develop even in the U' > U model''^. Therefore, if we 
try to explain the Cge softening within the RPA, we have 
to assume the case (ii), i.e., the cooperative Jahn- Teller 
transitionas, as proposed by Yanagi et alf^^ . However, it 
contradicts with experiments facts, that is, (a — 6)/(a + 
h) < 0.3% and Uxz — nyz S 5% in the orthorhombic 
phase irrespective of higher Ts (~ lOOK), and 6* > in 
the under-doped pnictides. Therefore, the RPA analysis 
done by Yanagi et al.'^^ is inconsistent with experiments. 

In this paper, we studied the structure transition in 
case (i), i.e., the quadrupole-quadrupole interaction type. 
The origin of quadrupole-quadrupole interaction is the 
"two-orbiton process" , which is not taken into account 
in the RPA. In two-dimensional systems, the two-orbiton 
process gives the Curie- Weiss behavior in Eq. (A2), if al- 
lowed by the orbital selectrion rule discussed in Sec. VI B. 
By using reasonable sets of parameters (^ac < g), experi- 
mentally observed Curie- Weiss behavior of Cee in under- 
doped pnictides {Ts > O)'^^'^'^ is well reproduced, given in 
Eq. (55): According to Eq. (54), Tafq (= S) is positive 
in under-doped pnictides, and Tafq ~ at the critical 
point Ts = 0, consistently with experimental reports"^^. 
In over-doped pnictides, the softening of Cge becomes 
moderate since the two-orbiton process gives very weak 
temperature dependence. 



Appendix B: Quadrupole susceptibilities in the 
ten-orbital model: effect of unfold gauge 
transformation 
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FIG. 14: (Color online) Quadrupole susceptibilities in the 
ten-orbital model for (a) Xxz{q), (b) X?z{q), and (c) XxyiQ), 
respectively. We put n — 6.05, T = 0.05, and Qc = 0.98. 








(n.O) 



FIG. 15: (Color online) Quadrupole susceptibilities in the ten- 
orbital model for (a) x^i^^iQ) and X^i^^il), (b) X?i^^{q) 
and Xyi^'^iQ), and (c) X?/^(<7) and X^y (q), respectively. 
We put n = 6.05, T = 0.05, and = 0.98. 



In Sec. Ill, we have calculated Xp(q,0) based on the 
five-orbital model using the RPA, shown in Fig. 4. In this 
Appendix, we discuss the effect of the unfold-gauge trans- 
formation in deriving the five-orbital model on Xp (q, 0). 
This gauge transition changes the signs of 2,3-orbitals 
for Fe-B sites. Therefore, signs of quadrupole opera- 
tors , at Fe-B sites are reversed, as recognized in 

xz/yz ' o 

Eqs. (4)-(8). For this reason, Xxz,yz{Qj^) is not gauge- 
invariant since it contains linear terms with respect to 
'^xz/yz' although the el-el interaction in Eqs. (14)-(16) 
is gauge invariant. 

From now on, we calculate the quadrupole susceptibil- 
ities in the ten-orbital model. We describe the orbitals 
of Fe-A (Fe-B) ions as 1 ~ 5 (6 ~ 10). In the RPA, the 
susceptibilities are given as 



^Q^AMAB,(^^) = E E °rX«'w('?)or'"', (B2) 



Q,ABi 



(Bl) 



//' mm' 



for r = 



xz,yz,xy, 



where Or 



o 



/ — 5,m — 5 



for Z, TO > 6. 



The obtained quadrupole susceptibilities in Eq. (Bl) is 
show in Fig. 14, and its diagonal and off-diagonal terms 
with respect to A and B are also shown in Fig. 15. In 
these figures, we have unfolded the susceptibilities into 
the single-iron BZ in order to make comparison with Fig. 
4. According to Eq. (Bl), Xr (?) ten-orbital model 
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is given by x?'^^{q) + Xr"^^(9)- We see that X^zi^^^) 
is Fig. 14 (a) has a sharp peak at q = (0,7r), since both 

Xxz'^^il) and Xxi'^^il) ill I^ig- 15 (a) have peaks at 
q= (0,7r). 

If we perform the unfold-gauge transformation, the 
sign of Xr'^^il) is inverted for F = xz and yz. There- 
fore, in the five-orbital model, Xril^^) — Xr'^^il) ~ 



Xp ' {q) for r = xz and yz. For this reason, the peak in 
Xxz(9) 0) move from q = (0, tt) to (tt, 0) under the gauge 
transformation, consistently with the result in Fig. 4 (a). 
In contrast to X^^/^^^ (</) , XxyiQ) is gauge invariant. For 
this reason. Fig. 14 (c) in the ten-orbital model coincides 
with Fig. 4 (c) in the five-orbital model. 
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